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Abstract 

The present study investigates linear and volatile (nonlinear) correlations of first- 
order autoregressive process with uncorrelated AR (f) and long-range correlated 
CAR (f ) Gaussian innovations as a function of the process parameter (9). In the light 
of recent findings [1], we discuss the choice of CAR (1) in modeling daily temperature 
records. We demonstrate that while CAR (1) is able to capture linear correlations it 
is unable to capture nonlinear (volatile) correlations in daily temperature records. 
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1 Introduction 



Noise at the dynamical and measurement level often encourages the choice of 
stochastic interpretation of time series data. Detrended fluctuation analysis 
(DFA) [2, 3] and its extensions [4-6] have been used widely to determine the 
nature of correlations from stochastic processes obtained from a wide-range of 
systems [7]. These studies have subsequently suggested the choice of stochastic 



* R.B. Govindan, Graduate Institute of Technology, University of Arkansas at Little 
Rock, Little Rock, AR 72204 

Email address: rbgovindan@git.ualr.edu (R.B. Govindan). 



Preprint submitted to Elsevier Science 



2 February 2008 



models to explain complex behavior observed in synthetic and natural data 
sets [1]. The objective of the present study is to understand the scaling be- 
havior of the first-order stochastic processes and their volatility series [4] as 
function of the process parameter. Based on the results obtained for the var- 
ious parameter regime, we evaluate the choice of these processes in modeling 
real- world temperature data [1]. 



2 First order auto-regressive process and its extension 

2.1 AR (1) with uncorrelated Gaussian innovations 

Classical AR (1) is given by the expression 

x n = 6>x n _i + e n ,n = l.JV, (1) 

where e n represents independent and identically distributed (i.i.d) Gaussian 
innovations with zero mean and unit variance. The corresponding volatility 
series of (1) is given by y n = \x n —x n -i\, where |iu| represents the absolute value 
of w. Correlations in y n imply that clusters of (small) big changes are likely 
to be followed by clusters of (small) big changes. It can be shown analytically 
that Eqn. (1) is stationary when \8\ < 1 i.e. — 1 < 8 < 1. The autocorrelation 
function (ACF) of (1) decays exponentially and is of the form p(k) = 8 h 
for lag k. The nature of the decay is governed by the process parameter 8. 
While p(k) = 8 k exhibits an exponential decay for (0 < 8 < 1), it exhibits an 
exponential decay with oscillations for (—1 < 8 < 0). It is important to note 
that (—1 < 8 < 0) and (0 < 8 < 1) represent cases where the current sample is 
negatively and positively correlated with its immediate past, respectively. The 
strength of the short-term correlated and short-term anti-correlated behavior 
is dictated by the process parameter 8. 



2.2 AR (1) with long-range correlated Gaussian innovations, CAR (1) 

We also consider the case where e n in (1) is generated by a long-range corre- 
lated Gaussian process [8]. This extended AR (1) with long-range correlated 
Gaussian innovations shall be referred to as CAR (1) in the subsequent sec- 
tions. The choice of the term CAR was encouraged by recent studies [1]. In 
the light of recent findings [1], we critically evaluate the choice of CAR (1) for 
modeling real-world daily temperature records. 
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3 Results and discussion 



A complete description of DFA and its extensions can be found elsewhere 
[3-5]. In the present study, the length of the data sets is chosen sufficiently 
larger N = 2 14 , in order to avoid bias in the estimation procedure due to finite 
sample-size effects. DFA with fourth order polynomial detrending is used in 
order to minimize bias due to local polynomial trends [4, 9]. We also integrate 
the data prior to DFA estimation and scale the resulting fluctuation function 
by the window size [4]. The fluctuation function of uncorrelated noise a = 0.5 
and the long-range correlated noise a = 0.8 and a = 0.65 shall be used as 
reference in the subsequent discussions. 



3.1 AR (1) with uncorrelated Gaussian innovations 



As noted earlier (Sec. 2) AR (1) process is stationary for — 1 < 9 < 1. In 
the present study, we consider parameters (9 = —0.1 and 0.1) and (9 = 
—0.93 and 0.93) . While the former is close to the uncorrelated regime (0 = 0), 
the latter is close to the non-stationary regimes, 9 = — 1 and 9—1. Log- log 
plots of the fluctuation function versus time scale for AR (1) and its volatil- 
ity series with 9 = 0.1 and — 0.1, Fig. la, is parallel to exponent 0.5 in the 
asymptotic regime, with a slight deviation from 0.5 in the short time scales for 
9 = —0.1. A similar analysis was carried for process 9 = 0.93 and — 0.93 Fig. 
lb. Unlike 9 = ±0.1 whose fluctuation function is relatively homogenous, 
characteristic crossovers are observed for 9 = ± 0.93 and their volatility series, 
Fig. lb. For 9 = 0.93, the log-log plot of the AR (1) process displayed cor- 
related behavior at short-time scales and uncorrelated behavior at the larger 
time scales a ~ 0.5. However, for 9 = —0.93 , one observes a characteristic 
crossover from anti-correlated behavior for time scales (s < 30) to uncorre- 
lated behavior for (s > 30). The volatility series for 9 = 0.93 was homogenous 
and exhibited an exponent of (a ~ 0.5) characteristic of monofractal data [4], 
whereas that of 9 = —0.93, displayed a marked crossover from correlated to 
uncorrelated behavior. 

From the above discussion, it is clear that the process parameter 9 has a signif- 
icance impact on the scaling behavior at shorter time scales. This is reflected 
by marked distortions and cross-overs in the fluctuation plots. It can also be 
noted that irrespective of the choice of 9, the scaling in the asymptotic regime 
resembles that of an uncorrelated noise, a = 0.5. 
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3.2 AR(1) with long-range correlated Gaussian innovations, CAR (1) 

Long-range correlated Gaussian innovations e n with scaling exponent of (a = 
0.8) [8] was used to generate CAR (1). The fluctuation function of CAR (1) 
process for 9 = ±0.1 were parallel to each other and to e n (a = 0.8), see Fig. 2a. 
The corresponding volatility series displayed an exponent of 0.5 characteristic 
of monofractal data. Thus the scaling behavior of CAR (1) is similar to that 
of AR (1) for process parameter 9 close to zero. 

For 9 = —0.93, at short time scale s < 35, CAR (1) is highly anti-correlated 
however the behavior mimicked that of e n (a = 0.8) for s > 40, Fig. 2b. For 9 = 
0.93, short time scales had an exponent a ~ 1.75 close to the integrated version 
of e n (a = 0.8). However, in the asymptotic regime the scaling exponent was 
that of e n (a = 0.8), Fig. 2b. Fluctuation function obtained for the volatility 
series (9 = 0.93) was parallel to the reference line with slope 0.5 characteristic 
of monofractal data. However, the fluctuation function of the volatility series 
(9 = —0.93) displayed correlated behavior at short time scales (s < 200) and 
uncorrelated behavior (a ~ 0.5) in the asymptotic regime. As in the case 
of AR (1), the process parameter 9 has a significant impact on the scaling 
behavior of CAR (1) process, revealed by marked distortion and crossovers in 
the log-log plots. 

3. 3 Critical note on the choice of CAR (1) to model temperature correlations 

It has been shown [10,11] that atmospheric temperature data (from a ran- 
domly chosen 14 continental regions) display power correlations with DFA 
exponent a = 0.65 and this is considered as universal persistence behavior 
of the temperature data. Using this result as a benchmark, the quality of the 
temperature data simulated by the Global Climate Models has been evaluated 
[12]. In a later study [13] it has been shown that temperature data from conti- 
nental and costal zones display lesser deviations from the previously observed 
universal behavior (exponent) of 0.65 while the islands display greater vari- 
ability with the average exponent being close to 0.8 [14]. In a recent study, [1] 
CAR(l) process was proposed as a plausible model for explaining correlations 
in atmospheric daily temperature fluctuations collected from fourteen meteo- 
rological centers in Hungary. The choice of the CAR (1) model was inspired 
by the fact that temperature records exhibit both short-term and long-term 
correlations [1, 7]. We generated the data using the CAR (1) model proposed 
by in [1] 

Xi = («i - c)xi-i + 5rji. (2) 

The model parameters were chosen as a\ = 0.8, 5 = 2.1, r\i is zero-mean, unit 
variance long-range correlated Gaussian innovations with scaling exponent 
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a = 0.65, p = a — 0.5, c = 2p 3 / 2 as in [1]. 



Log- log plot of the fluctuation function of Eqn. 2, Fig. 3, displayed strong 
correlations with exponent (a = 0.85) at short scales (s < 100). However, 
in the asymptotic regime it displayed an exponent (a ~ 0.65) close to that 
of the long-range correlated Gaussian innovations. Fluctuation function of 
the corresponding volatility series was parallel to (a = 0.5), characteristic 
of monofractal processes. This has to be contrasted with recent studies [7], 
which provided compelling evidence of long-range correlations (a = 0.6) in 
the volatility series of daily temperature fluctuations. Such long-range volatile 
correlations were observed across randomly chosen and spatially separated me- 
teorological sites across the world. This in turn was attributed to universality 
in daily temperature records [7]. 

It has been shown that long range correlated data with exponents 0.5 < a < 
1.5 will exhibit uncorrelated behavior in volatility analysis while the multifrac- 
tal signals exhibit long range correlations in volatility analysis [4, 6] . As the 
CAR (1) model, Eqn. 2, is driven by a linear process (see above), one would ex- 
pect the model to display uncorrelated behavior in the volatility analysis and 
is indeed the case obtained for this model. Since, the CAR (1) model Eqn. 2 
[1] captures the linear features of the temperature data (as the volatility ex- 
ponent is 0.5 for this model), it will be unable to account for the nonlinear 
features in the data irrespective of the choice of the process parameter, as 
demonstrated by Figs. 2 and 3. Equation (2) [1] does show distinct scaling 
behavior with a marked crossover, however, it is not appropriate for modeling 
the dynamics of daily temperature records. In [1], the authors subsequently 
proposed a nonlinear counterpart of CAR(l) model, termed NLCAR(l) in or- 
der to capture the probability distribution of the daily temperature records. 
The DFA results obtained in [1] for NLCAR(l) and CAR(l) were indistin- 
guishable (for details we refer to [1]). However, this does not necessarily imply 
NLCAR(l) and CAR(l) have similar volatile correlations. Recent studies [15] 
have used nonlinear stochastic models successfully to model glacial dynamics 
also reflected in long-range volatile correlations. Encouraged by previous find- 
ings [4,7,15] and results presented in [1], we believe that NLCAR(l) may be 
a more appropriate model to capture dynamics of daily temperature records. 
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Fig. 1. Scaling analysis of AR (1) process, (a) Fluctuation function for 9 = —0.1 
and 9 = 0.1 is shown by (solid curve with open circles) and (solid curve with 
x) respectively. Corresponding volatility analysis for 9 = —0.1 and 9 = —0.1 are 
shown with (open circles and x) respectively, (b) Fluctuation function for 9 = —0.93 
and 9 = 0.93 is shown by (solid curve with open circles) and (solid curve with x) 
respectively. Corresponding volatility analysis for 9 = —0.93 and 9 = 0.93 are shown 
with (open circles and x) respectively. Fluctuation function of the uncorrelated 
noise with slope 0.5 is shown at the bottom of Fig 1 (a) and (b) as reference. In 
the asymptotic regime all the curves are parallel to the reference line indicating 
uncorrelated nature of the process. 
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Fig. 2. Scaling analysis of CAR(l) process (i.e. AR(1) process driven by long range 
correlated noise with DFA exponent 0.8) for process parameters (9 = ±0.1 and 
= ±0.93) and their volatility series is shown in (a) and (b) respectively. Expla- 
nations are same as in Fig. 1. Fluctuation function of the uncorrelated noise (slope 
0.5, solid line) and long-range correlated noise (slope = 0.8, dotted line) is shown 
at the bottom in (a) and (b) as reference. In the asymptotic regime DFA curves 
for the original data are parallel to the reference line with slope 0.8 and the DFA 
curves of the volatile data are parallel to the reference to line with slope 0.5. 
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Fig. 3. Scaling analysis of CAR(l) temperature model (Sec. 3, Eqn. 2). The fluctu- 
ation function for the temperature model is shown in solid curve with open circles. 
The fluctuation function of the corresponding volatility series is shown with open 
circle. Fluctuation function of the uncorrelated noise (slope 0.5, solid line) and 
long-range correlated noise (slope = 0.65, dotted line) is shown at the bottom in 
(a) and (b) as reference. In the asymptotic regime, DFA curve of the original data 
is parallel to the reference line with slope 0.65 while at short time scales s < 200 
the data are correlated with an exponent a = 0.85. The DFA curves of the volatile 
data are parallel to the reference line with slope 0.5 indicating uncorrelated nature 
of the volatile correlations. 
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